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Abstract. We report on a new force that acts on cavities (literally empty regions of 
space) when they are immersed in a background of non-interacting fermionic matter 
fields. The interaction follows from the obstructions to the (quantum mechanical) 
motions of the fermions caused by the presence of bubbles or other (heavy) particles 
in the Fermi sea, as, for example, nuclei in the neutron sea in the inner crust of a 
neutron star or superfluid grains in a normal Fermi liquid. The effect resembles the 
traditional Casimir interaction between metallic mirrors in the vacuum. However, the 
fluctuating electromagnetic fields are replaced by fermionic matter fields. We show that 
the fermionic Casimir problem for a system of spherical cavities can be solved exactly, 
since the calculation can be mapped onto a quantum mechanical billiard problem of a 
point-particle scattered off a finite number of non-overlapping spheres or disks. Finally 
we generalize the map method to other Casimir systems, especially to the case of a 
fluctuating scalar field between two spheres or a sphere and a plate under Dirichlet 
boundary conditions. 



PACS numbers: 03.65.Nk, 03.65. Sq, 03.75.-b, 21.10.Ma, 74.45. +c 
1. Introduction 



1.1. The original Casimir effect 

In 1948 the Dutch physicist H.B.G. Casimir predicted the remarkable effect 1] that two 
parallel, very closely spaced, uncharged metallic plates attract each other in vacuum. 
The origin of this force can be traced back to the altered (mode sums over the) zero- 
point fluctuations of the electromagnetic field which are induced by the presence of the 
two plates which are added relative to the free case, or rather, which are brought from 
very far separation into a close distance. The distinctive property of the Casimir effect 
is that its strength and, perhaps, its sign are geometry-dependent (for a review see |2j). 
Our aim is to generalize two features of the Casimir effect, the mode sum structure 
and the geometry-dependence, to the fermionic Casimir effect, where the fluctuating 
photons are replaced by (non-relativistic) Fermi fields. 
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1.2. Utilizing the geometry dependence of the Casimir energy 

Let us invert the logic and define the Casimir energy as the energy resulting from the 
geometry- dependent part of the density of states (d.o.s.) - a concept that is closely 
related to the shell correction energy in nuclear physics: 

p(E) = E k ) = po(E) + Pbulk(£) + 5p C (E, geom.-dep.) . (1) 

Here {Ek} are the eigenenergies of the modes, po is the d.o.s. of the homogeneous 
background and Pbuik is the bulk d.o.s. that sums up the excluded volume effects, surface 
contributions and Friedel oscillations caused by each of the obstacles separately. The 
remaining part Spc is of central interest to us. It is the only term which "knows" about 
the relative geometry-dependence of the obstacles. Now the Casimir energy can be 
extracted from the geometry- dependent part of the density of states as a simple integral 

Sq = J dE E 5 pc(E, geom.-dep.) = — J dE Afc(E, geom.-dep.) , (2) 

where Mc(E) is the geometry-dependent part of the integrated density of states or 
number of states, N(E) = Ee, ®(E - E k ) = / E dE' p(E'). 

1.3. Generalization of the Casimir energy concept to matter fields 

Let us assume that space is not "filled" with fluctuating electromagnetic modes, but 
with a gas of non-interacting (non-relativistic) fermions. Under this scenario we have the 
following similarities with the ordinary Casimir effect: In both cases, there exist mode 
sums, J2k huJki with constant degeneracy factors. The constant weight of the fluctuating 
electromagnetic modes can be traced back to the two helicity states of the photon. The 
constancy in the degeneracy of the fermionic matter modes, on the other hand, follows 
from Pauli's exclusion principle, where the pertinent weight can be formulated in terms 
of spin and isospin factors. However, the Casimir mode summation over matter fields 
differs from the case of fluctuating fields by the presence of a new independent scale in 
addition to the geometric scales (e.g. the separation L and area A of the plates), namely 
by the presence of the Fermi energy (i.e., the chemical potential p at zero temperature). 

In the following we will consider the case of matter fields (non-relativistic fermions) 
located in the space between voids or cavities, such that the matter fields will build up 
a quantum pressure on the voids. Even if we assume that the matter fields are non- 
interacting and the voids are non-overlapping, an effective interaction between these 
empty regions of space will still arise in the background of the non-interacting fermionic 
matter fields, since the cavities - depending on their geometric arrangement - can shield 
the free movement of the matter modes such that a net change in the mode sum over 
the Fermi states in the Fermi sea results. 

Applications of this scenario exist e.g. for the inner crust of neutron star. With 
increasing distance from the star's surface (inside the star), the nuclei start to loose 
neutrons due to the growing pressure and density. Bulgac and Magierski argued that 
the shell correction energies (in other words the fermionic Casimir energies) of the 
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resulting low-density neutron matter in the presence of bubbles (i.e. nuclei) are of the 
same order of magnitude [3] as the differences between the usually assumed liquid-drop- 
model sequence of phases jl] . Thus disordered lattices of bubbles are in competition with 
the standard sequence of (i) nuclear drops (nuclei), (ii) nuclear rods, (hi) nuclear plates 
(all three surrounded by a low-density background of neutrons), (iv) tubes, (v) bubbles 
(both surrounded by the high-density background of nuclear matter), until finally a 
uniform nuclear matter phase is reached. 

An analogous investigation of bubbles inside a Fermi gas background is also of 
relevance for the inner core of neutron stars. Under the assumption that quark matter 
does exist and quark droplets can form there, a similar pattern is predicted, with quark- 
droplets (bags) inside hadronic matter taking over the role of the embedded nuclei. 

Also in the laboratory, the study of the interaction of cavities inside a uniform 
fermionic background could be of importance Examples are Ceo buckyballs immersed 
in liquid mercury. The liquid metal itself serves only as a non-rigid neutral background 
which provides the Fermi gas environment via its conductance electrons, in which the 
buckyballs "drill" the voids. Another example would be buckyballs in liquid 3 He as 
Fermi gas. Finally, in the future, boson-condensate droplets immersed in dilute atomic 
Fermi condensates could serve as systems with which the effective interactions of cavities 
inside a Fermi gas could be studied in the lab. 

2. Casimir calculation mapped onto a scattering problem 

Note that the Casimir calculation for fermionic (non-relativistic) matter fields simplifies 
enormously, because the presence of a second scale, the chemical potential p=h 2 k 2 F /2m 
(or the Fermi momentum kp), provides for a natural UV-cutoS, Auv = A* (here m is the 
mass of the fermion). Thus the Casimir energy for fermions between two impenetrable 
(parallel) planes at a distance L is simply given by the chemical potential times a finite 
function of the dimensionless argument kpL, £q = pF(kpL). 

For more complicated geometries, the computations become more and more 
involved as it is the case for the ordinary electromagnetic Casimir effect. However, 
Casimir calculations of a finite number of immersed non-overlapping spherical voids or 
rods, i.e. spheres and cylinders in 3 dimensions or disks in 2 dimensions, are still doable. 
In fact, these calculations simplify because of Krein's trace formula |EIIZ| 

Sp(E) = p{E) - p (E) = -L ^-tvlnS n (E) , (3) 

2m ahi 

which links the variation in the level density Sp(E) (the difference of the total density of 
states and the background one) to the energy- variation of the phase shift ~kid.e&S n (E) 
of the n-sp here/disk scattering matrix S n (E). Note that the level densities on the left 
hand side are averaged over an energy-interval larger than the mean-level spacing in the 
volume V of the entire system in order to match the continuous expression on the right 
hand side. In this way the Casimir calculation is mapped to the quantum mechanical 
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analog of a classical "billiard" problem: the hyperbolic or even chaotic scattering of a 
point-particle off an assemble of n non-overlapping spheres (or disks) [HI El HII1HH H2 UH]- 
Moreover, it is now possible to extract the geometry-dependent Casimir fluctuations 
from the multiple-scattering part of the scattering matrix. In fact, as shown in 
[TT| H21 US], the determinant of the n-sphere/disk S-matrix can be separated into a 
product of the determinants of the 1-sphere/disk S-matrices S\(E,cii), where a, is the 
radius of the single scatterer i, and the ratio of the determinant of the inverse multi- 
scattering matrix M(k) and its complex conjugate: 

, n f-A-i rwr, , ] det (M(fc*)) f 

detS n (E) = jndetftCfl.aOj d( ^(fc) ■ (4) 

For the fermionic Casimir case, the energy E and the wave number k are related as 
E = h 2 k 2 /2m. When inserted into Krein's formula, the product over the single-scatterer 
determinants generates just the bulk (or Weyl term) contribution to the density of states 

n 1 d n 

Pbulk(£, H}) = E/°Weyl(£,ai) = -^—Jfi E lndet5 'l(' £ '> a i) ' ( 5 ) 
2=1 i=l 

which takes care of the excluded volume terms and the surface terms (including Friedel 
oscillations). The geometry-dependent part of the d.o.s. is therefore given by a modified 
Krein equation [3] which is formulated in terms of the (inverse) multi-scattering matrix 
instead of the full S-matrix 

n 

5p c (E, {ai}, {r i:j }) = p(E) - p (E) - E P Wcy i(£, a<) 

i=l 

= -^Im |^lndetM[fc(E),{a,},{rl,}]| , (6) 

where r^- are the relative separation vectors between the centers of the spheres (or disks). 
The pertinent Casimir energy can then be read off from the finite integral 

8 C = r dE (E - fi) 5p c = - dE Jf c . (7) 
Jo Jo 

3. The calculation 

Equation (jB)) does allow to simplify the problem, since there exists a close-form 
expression for the inverse multi-scattering matrix for n spheres (of radii dj and mutual 
separation r^y, labelled by the indices j,j' = 1, 2, • • ■ , n) in terms of spherical Bessel and 
Hankel functions, spherical harmonics and 3j-symbols, where I, I' and m, m' are total 
angular momentum and pertinent magnetic quantum numbers, respectively j!3j : 



MH iVm , = 5"'5 w 5 mml + (1 - i 2m+l '- 1 v /4tt(2/ + 1)(2/'+1) 

WJ hP(ka f ) kom^-v V ) 

x ( J. m n J l m )^UJ^{kr jr )Y^'(f^). (8) 
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The unit vectors r^-> point from the origin of sphere j (as measured in its local coordinate 
system) to the origin of sphere j'. The local coordinate system of sphere j' is mapped 
to the one of sphere j with the help of the rotation matrix D l mlih {j^j'). 

For small scatterers the expression of the multi-scattering matrix even simplifies: 

M jj '(E) » 5 jj ' - (1 - <F')//(fl) ^ + OOwave). (9) 

since only s-wave scattering is important; i.e., spherical waves modulated by s-wave 
amplitudes fij(E) propagate between the spheres. The integrated d.o.s. in the case of 
two small spherical cavities of common radius a and center-to-center separation r is [H] 

i n 2 

Af°°(E) = ImlndetM 00 (£) « z/ deg — sin[2(r - a)k] + O Ukaf) (10) 

7r irr 2 ^ ' 

where ^d eg is the spin/isospin-degeneracy factor. This expression should be compared 

with the semiclassical approximation that sums up all partial waves 

M° C U E ) = v *° H-Kr{r-2a) ~ 2a) * ] ' (U) 

In fact, the latter is the leading contribution to Gutzwiller's trace formula [Hj, namely 
the non-repeated contribution of the two-bounce periodic orbit between the two spheres, 
with the action S po (k)/h = 2(r — 2a)k where 2(r — 2a) is the length of the geometric 
path. Note that the semiclassical result is suppressed by a factor of 1/4 in comparison 
to the small-scatterer one. 

As shown in Ref. |Hj the semiclassical result is a very good approximation of the full 
quantum mechanical result calculated from the exact expression (JHJ) of the two-sphere 
scattering matrix when plugged into the modified Krein formula (jOJ. Therefore the 
Casimir energy for the two spherical cavities inside a non-relativistic non-interacting 
fermion background can be approximated in terms of a spherical Bessel function j\ as 

£°° = dE^S°(E) » -^ deg // Q ji[2(r - 2a)k F ] , (12) 

which is valid for kpa > 1. This expression is long-ranged, i.e. 1/L 3 with L = r — 2a. 
For the sphere-plate system the Casimir energy reads instead 

£c ~ -^dcg /^ttt" \3ii 2 ( r ~ a ) k F] , (13) 

2ir(r — a) 

which scales even as 1/L 2 with L = r — a. Note that in both cases, the two-sphere system 
or the sphere-plate system, the Casimir energy does not have a fixed sign in contrast 
to the standard Casimir effect with fluctuating electromagnetic or scalar fields between 
these obstacles. Instead the sign of the Casimir energy oscillates as function of the action 
of the two-bounce orbit. Therefore, under increase of the distance between the cavities, 
the Casimir energy, which starts out to be attractive, can be made repulsive, and under 
a further increase of this distance, it can become attractive again, with a decreased 
strength of course. The reason for this new type of behaviour of a Casimir energy is the 
presence of a new scale in addition to the length scales, namely the chemical potential 
fi. In fact, the strength of this fermionic Casimir energy scales with the strength of the 
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chemical potential and therefore with the UV-cutoff of the theory. Also this behavior 
distinguishes the fermionic Casimir effect from the standard Casimir effect: the latter 
is governed by the infrared behavior of the corresponding density of states. 

4. Fermionic Casimir effect between superfluid grains 

A generalization of the fermionic Casimir effect to the case of superfluid cavities 
immersed in normal fermionic matter is reported in Ref. ^2]- The pairing interaction 
due to the superfluid obstacles leads to enhanced Casimir contributions, because 
of the dominance of the particle-hole terms over the particle-particle and hole-hole 
contributions. Semiclassically, this can be explained by the focusing (and only for large 
separations defocusing) nature of the Andreev reflections fB| (see Fig. [TJ in comparison 
to the defocusing specular reflections at normal circular boundaries. Since the particle- 




N 



Figure 1. Si and S2 are two superfluid grains immersed in a normal Fermi liquid N. 
Indicated is the Andreev reflection between an incoming particle (hole) and an outgoing 
hole (particle). The figure is from Ref. |15| . 

hole d.o.s has the opposite sign to the particle-particle one and since the d.o.s has to 
be integrated over the quasi-particle energy as measured from the Fermi surface, the 
resulting Casimir energy for superfluid obstacles is strongly repulsive [TK] . 

5. Application of the map method to the scalar Dirichlet problem 

Note that the map method onto a scattering problem can be generalized from the 
Casimir effect in the Fermi sea to other systems; it is especially applicable for the case 
of a fluctuating scalar field between two spheres or a sphere and a plate (with Dirichlet 
boundary conditions) [T71 ll8|. Whereas the s-wave approximation in the small-scatterer 
limit is not important for the fermionic Casimir energy which is governed by the UV 
part of the d.o.s. (i.e. by the contribution at the chemical potential), it is essential for the 
large distance physics of the fluctuating-scalar Casimir effect, as the latter is governed by 
the infrared behaviour of the d.o.s. Remember the relative factor of four between the s- 
wave and semiclassical result of the two-cavity d.o.s. The same expressions are valid for 
the d.o.s. of the fluctuating-scalar Casimir effect. Therefore, the Casimir energy for two 
far separated Dirichlet spheres is enhanced by a factor of 4 x (90/-7T 4 ) ^H] relative to the 
result of the proximity force approximation which is only valid for small separations [20J . 
The extra 90/7T 4 takes into account the missing repeats of the two-bounce orbit in the 
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semiclassical formula which are unimportant at large separations, but needed at 
small ones (see curve D of Fig- EJ) - For the sphere-plate configuration the corresponding 
enhancement at large separations is 2 x 90/7T 4 , because only one sphere is involved; see 
Fig.Elfor the comparison of this result (C) with the exact calculation (A) and the s-wave 
approximation (B) jTHI and with the proximity force approximations of Refs. [T7] I18j. 
The exact data (A) are compatible with numerical data of the wordline approach, within 
error bars, ^7] in the window of overlap, i.e. L < 4a. In the case of the electromagnetic 
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Figure 2. Predictions for the scalar Casimir energy £c{a,L) of the Dirichlet sphere- 
plate configuration in units of — IAOOTictt 3 a / L 2 as function of the distance L divided by 
the radius a of the sphere. The circles (A) represent the numerically calculated exact 
expression for the sphere-plate system, the curve (B) shows the s-wave approximation, 
the line (C) represents the asymptotic limit 2 x 90/ir 4 , the crosses (D) represent the 
numerically calculated results from the semiclassical Gutzwiller formula including all 
repeats, see [T^J- The curves (E) denotes the plate-based, the curve (F) the sphere- 
based proximity- force approximation, see |17l 118) . 

Casimir effect, the s-wave dominance at large separation has to be replaced by a p-wave 
dominance, since the charge-neutrality of the sphere forbids a monopole term, whereas 
the standard Casimir-Polder energy is dominated by induced-dipole contributions. 



6. Summary 

We have shown that there exists an effective interaction between voids inside a Fermi gas 
background, even if the fermions are non-interacting. This new form of Casimir energy 
is neither attractive nor repulsive, but oscillates according to the relative arrangement of 
the cavities (in the superfluid case it is strongly repulsive). The spheres can be replaced 
with other objects, if the curvature radii are larger than the Fermi wave length. The 
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effects of finite surface thickness can be booked as Weyl-term contributions and do not 
affect the geometrical Casimir part if the objects do not overlap. Since the Casimir 
interaction between bubbles is oscillating and rather long-ranged (in comparison to van 
der Waals terms), disordered lattices are expected as emerging structures. The disorder 
can be further enhanced by finite temperature T and corrugated surfaces. The map 
method can be generalized to other Casimir systems, especially to the scalar Casimir 
effect for Dirichlet spheres or sphere-plate systems. 
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